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Abstract 

The angular distributions of the baryon-antibaryon low-mass enhancements seen in 
the charmless three-body baryonic B decays B + —> ppK + , B° —> ppKg, and B° —> 
pAtt~ are reported. A quark fragmentation interpretation is supported, while the 
gluonic resonance picture is disfavored. Searches for the 0 + and 0 ++ pentaquarks 
in the relevant decay modes and possible glueball states Q with 2.2 GeV/c 2 < 
M p p < 2.4 GeV/c 2 in the pp systems give null results. We set upper limits on the 
products of branching fractions, B(B° —> Q + p) x B(Q + —» pKg) < 2.3 x 10” 7 , 
B(B + -► 0++p) x B(Q ++ -»• P K +) < 9.1 x 10“ 8 , and B(B + -»• QK+) x B(G -»• 
pp) < 4.1 x 10”" at the 90% confidence level. The analysis is based on a 140 fb” 1 
data sample recorded on the T(4S) resonance with the Belle detector at the KEKB 
asymmetric-energy e + e” collider. 

PACS: 13.25.Hw, 13.60.Rj 
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Observations of several baryonic three-body B decays have been reported re¬ 
cently [1,2,3,4], One common feature of these observations is the peaking of 
the baryon-antibaryon pair mass spectra toward threshold, as originally con¬ 
jectured in Refs. [5,6] and elaborated more recently in Refs. [7,8,9]. The same 
peaking behavior near threshold has been found in baryonic J/ij) decays [10] 
as well, indicating that this may be a universal phenomenon. Possible ex¬ 
planations include intermediate (gluonic) resonant states or non-perturbative 
QCD effects of the quark fragmentation process [7,8]. Alternatively, the dy¬ 
namical picture can be replaced by an effective range analysis with a baryon 
form factor [9]. To distinguish among the above hypotheses for the production 
mechanism, we study the angular distributions of the threshold enhancements 
in the helicity frame for the B + —> ppK + , B° —> ppKg and B° —> pA.n~ [11] 
modes. 

We use a 140 fb _1 data sample, consisting of 152 xlO 6 BB pairs, collected 
with the Belle detector at the KEKB asymmetric energy e + e _ (3.5 on 8 GeV) 
collider [12]. The Belle detector is a large solid angle magnetic spectrometer 
that consists of a three layer silicon vertex detector (SVD), a 50 layer central 
drift chamber (GDC), an array of aerogel threshold Cerenkov counters (ACC), 
a barrel-like arrangement of time of flight scintillation counters (TOF), and 
an electromagnetic calorimeter comprised of CsI(Tl) crystals located inside a 
superconducting solenoid coil that provides a 1.5 T magnetic field. An iron 
flux return located outside of the coil is instrumented to detect K° L mesons 
and to identify muons. The detector is described in detail elsewhere [13]. 

The event selection criteria are based on the information obtained from the 
tracking system (SVD and CDC) and the hadron identification system (CDC, 
ACC, and TOF). All primary charged tracks are required to satisfy track 
quality criteria based on the track impact parameters relative to the interac¬ 
tion point (IP). The deviations from the IP position are required to be within 
±1 cm in the transverse (x-y) plane, and within ±3 cm in the z direction, 
where the z axis is opposite the positron beam direction. For each track, the 
likelihood values L p , Lk, and L n that it is a proton, kaon, or pion, respec¬ 
tively, are determined from the information provided by the hadron identi¬ 
fication system. The track is identified as a proton if L p / (L p + Lk) > 0.6 
and L p /(L p + L n ) > 0.6, or as a kaon if L K /{L K + L n ) > 0.6, or as a pion 
if L n /(L K + L n ) > 0.6. The proton selection efficiency is about 84% (88% 
for p and 80% for p) for particles with momenta at 2 GeV/c, and the fake 
rate is about 10% for kaons and 3% for pions. Candidate Kg mesons are re¬ 
constructed from pairs of oppositely charged tracks (both treated as pions) 
having a mass consistent with the Kg nominal mass, \M 7T + 7T - — M K o| < 30 
MeV/c 2 , as well as a displaced vertex and flight direction consistent with an 
origin at the IP. Candidate A baryons are reconstructed from pairs of oppo- 
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sitely charged tracks—treated as a proton and negative pion—whose mass is 
consistent with the nominal A baryon mass, 1.111 GeV/c 2 < M pn - < 1.121 
GeV/c 2 . The proton-like daughter is required to satisfy L p /(L p + L n ) > 0.6. 

Candidate B mesons are reconstructed in the B + —> ppK + , B° —» ppKg, 
and B° —> pAn~ modes. We use two kinematic variables in the center of 
mass (CM) frame to identify the reconstructed B meson candidates: the beam 
energy constrained mass M bc = \J E^ eaia — p 2 B , and the energy difference A E = 
E B —E heam , where E he am is the beam energy, and p B and E B are the momentum 
and energy, respectively, of the reconstructed B meson. The candidate region 
is defined as 5.20 GeV/c 2 < M bc < 5.29 GeV/c 2 and —0.1 GeV < A E < 0.2 
GeV. From a GEANT [14] based Monte Carlo (MG) simulation, the signal 
peaks in the subregion 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 and \AE\ < 0.05 
GeV. The lower bound of A E is chosen to exclude possible contamination 
from so-called “cross-feed” baryonic B decays. 

The background in the fit region arises solely from the continuum e + e^ —> qq 
(q = w, d , s, c ) process. We suppress the jet-like continuum background 
events relative to the more spherical BB signal events using a Fisher discrim¬ 
inant [15] that combines seven event shape variables, as described in Ref. [16]. 
Probability density functions (PDFs) for the Fisher discriminant and the co¬ 
sine of the angle between the B flight direction and the beam direction in the 
T(4S) rest frame are combined to form the signal (background) likelihood C s 
(£{,). The signal PDFs are determined using signal MC simulation; the back¬ 
ground PDFs are obtained from the side-band data with M bc < 5.26 GeV/c 2 . 
We require the likelihood ratio 1Z = C s / ( C s + £ 5 ) to be greater than 0.7, 0.75, 
and 0.8 for the ppK + , ppKg, and pAtt~ modes, respectively. These selection 
criteria are determined by optimization of n s /^/n s + n b , where n s and n b de¬ 
note the expected numbers of signal and background events, respectively. We 
use the branching fractions from our previous measurements [2,3] in the cal¬ 
culation of n s . If there are multiple B candidates in a single event, we select 
the one with the best y 2 value from the vertex fit. 

We perform an unbinned likelihood fit that maximizes the likelihood function, 

p -(N s +N b ) N 

L = --- 1 ] [NsPsiM^AEi) + N b P b (M hCi , AE,)} , 

to estimate the signal yield in 5.20 GeV/c 2 < M bc < 5.29 GeV/c 2 and —0.1 
GeV < A E < 0.2 GeV; here P s ( P b ) denotes the signal (background) PDF, N 
is the number of events in the fit, and N s and N b are fit parameters representing 
the number of signal and background events, respectively. 

For the signal PDF, we use the product of a Gaussian in M bc and a double 
Gaussian in A E. We fix the parameters of these functions to values determined 
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by MC simulation [17]. The continuum background PDF is taken as the 
product of shapes in M bc and A E, which are assumed to be uncorrelated. 
We use the parameterization first used by the ARGUS collaboration [18], 
/(Mb c ) oc M bt yi — x 2 exp[—£(1 — a; 2 )], to model the M bc background, with x 
given by M bc /77 beam and / as a fit parameter. The A E background shape is 
modeled by a linear function whose slope is a fit parameter. 

The differential branching fraction as a function of the baryon pair mass 
is shown in Fig. 1. The error bars include the statistical uncertainty from 
the fit and the systematic uncertainty. Here, the efficiency as a function of 
baryon pair mass for each signal mode is determined by MC simulation. 
We sum these partial branching fractions to obtain: B(B + —» ppK + ) = 
(5.30lg;^ ± 0.58) x 10~ 6 , B(B° -> ppK° s ) = (1.20±g;ll ± 0.14) x 10~ 6 , and 
B(B° —» pAtt~) = (3.27/°;5i ± 0.39) x 10 -6 which are in good agreement with 
previous measurements [1,2,3]. These results also supersede our previous mea¬ 
surements with better accuracy. Note that we have imposed a charmonium 
veto for the ppK + and ppKg modes: the regions 2.850 GeV/c 2 < M pP < 3.128 
GeV/c 2 and 3.315 GeV/c 2 < M pp < 3.735 GeV/c 2 are excluded to remove 
background from B decay modes containing an r/ c , J/^>, //, y c0 , or Xci meson. 
The width of the low mass enhancement in each distribution of Fig. 1 depends 
on the signal mode. A narrow width is also observed in the newly discovered 
B + —> AA/l + decay [4], 

Systematic uncertainties are determined using high statistics control data sam¬ 
ples. For proton identification, we use a A —>• pn~ sample, while for K/tt 
identification we use a D* + —> D°n + , D° —> K~n + sample. Tracking effi¬ 
ciency is measured with fully and partially reconstructed D* samples. The 
uncertainty of Kg reconstruction due to off-lP tracks is determined from a 
D~ —> Kgn~ sample. The A and Kg reconstruction efficiencies have the same 
uncertainty due to off-IP tracks if the uncertainty of the daughter proton iden¬ 
tification criterion is not taken into account. The 7 Z continuum suppression 
uncertainty is estimated from control samples with similar final states, for 
example, B + —> J//>A' + with J/ijj —> p + /i^. Based on these studies, we as¬ 
sign a 1% error for each track, 3% for each proton identification, 2% for each 
kaon/pion identification, 5% for Kg and A off-IP reconstruction and 6% for 
the 7 Z selection. 

A systematic uncertainty of 4% in the fit yield is determined by varying the 
parameters of the signal and background PDFs. The MC statistical uncer¬ 
tainty contributes a 2% error in the branching fraction determination. The 
error on the number of BB pairs is 0.5%, where we assume that the branching 
fractions of T(4S) to neutral and charged BB pairs are equal. 

We first sum the correlated errors linearly and then combine them with the 
uncorrelated ones in quadrature. The total systematic uncertainties are 11%, 
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12%, and 12% for the ppK + , ppK° s , and pAn modes, respectively. 

We require the mass of the baryon pair to be less than 2.85 GeV/c 2 for 
the study of the threshold enhancement effect. The M bc distributions (with 
\AE\ < 0.05 GeV), and the A E distributions (with M bc > 5.27 GeV/c 2 ) for 
the ppK + , ppKg and pAn~ modes are shown in Fig. 2. The projections of the 
fit results are shown in Fig. 2 by solid curves. The B yields are 217 ± 17, 28.6 
and 48.8 l®; 2 for the ppK + , ppKg, and pAn~ modes, respectively. The 
measured branching fractions by summing the partial branching fractions in 
mass bins below 2.85 GeV/c 2 are B(B + —> ppK + ) = (4.59+o;g4 ± 0.50) x 10 -6 , 
B(B° -> ppK° s ) = (1.04ioli9±0.12) x 10~ 6 , and B{B° -»■ pAvr") = (2.62±g;^± 
0.31) x 10" 6 . 

We study the proton angular distribution of the baryon-antibaryon pair sys¬ 
tem in its helicity frame. The angle 9 P is defined as the angle between the 
proton direction and the meson direction in the baryon-antibaryon pair rest 
frame. Note that the angle is determined by p and K + (or p and K~) in the 
ppK + mode for definiteness. Fig. 3(a)-(c) shows the branching fractions as a 
function of cos 6% We define the angular asymmetry as A = where 

N + and N_ stand for the efficiency corrected B yield with cos 9 P > 0 and 
cos 9 P < 0, respectively. The angular asymmetry is determined to be 0.59jlolof 
for the ppK + mode. The asymmetry of the distribution indicates that the 
fragmentation picture is favored. Antiprotons are emitted along the K + di¬ 
rection most of the time, which can be explained by a parent b —> s penguin 
transition followed by su fragmentation into the final state as shown in Fig. 4. 
The energetic s quark picks up the u quark from a uu pair in vacuum and 
the remaining u quark then drags a ud diquark out of vacuum. This simple 
picture can describe the p — K + angular correlation. The spectator u quark 
and leftover ud diquark form an proton. 

The cos 9 P distribution of the ppK° s mode can not support nor refute this frag¬ 
mentation interpretation because of low statistics and no flavor information. 
The distribution for the pAn~ mode is quite flat (i.e. in favor of the gluonic pic¬ 
ture), in contrast to that of the ppK + mode, although both presumably share 
a common origin in the b —» s transition. In fact, this parentage suggests that 
it would be useful to examine the proton angular distribution relative to the A 
direction in the pn~ rest frame. We remove the M p \ < 2.85 GeV/c 2 constraint 
in order to check the full angular region; the result is shown in Fig. 3(d). It is 
evident that the fragmentation interpretation is supported: the proton tends 
to emerge parallel to the A baryon. 

As a cross check, the distribution of cos 9 P for background events in the ppK + 
sample is shown in Fig. 3(e). Similar distributions are obtained for the back¬ 
grounds of the ppKs and pAti~ modes. The background has a 1 + a cos 2 9 P 
distribution, which can be explained as arising from the random combina- 
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tion of two high momentum particles from the qq jets. The continuum MC 
simulation also confirms this feature. 

The newly observed narrow pentaquark state, © + [19], can decay into pK° s . 
We perform a search in B° —> ppK$ by requiring 1.53 GeV/c 2 < M pK o < 1.55 
GeV/c 2 . The M^ c and A E projection plots in Fig. 5 show no evidence for a 
pentaquark signal. Since there are few events in the fit window, we fix the 
background shapes from sideband data. We use the fit results to estimate the 
expected background (0.42 ± 0.13) and compare this with the observed one 
event in the signal region to set an upper limit on the signal yield of 3.9 events 
at the 90% confidence level [20,21], The systematic uncertainty is included in 
this limit. The related upper limit on the product of branching fractions is 
B(B° —> Q + p) x £>(© + —> pKg) < 2.3 x 10 -7 . We also perform a search 
for @ ++ , which can decay to pK + in the mode B + —> ppK + [22], Because 
there are only theoretical conjectures for the existence of such a state, we 
examine the wider mass region of 1.6 GeV/c 2 < M pK + <1.8 GeV/c 2 . We find 
no evidence for signal. Assuming this state is narrow and centered near 1.71 
GeV/c 2 , the upper limit on the yield is 3.3 events at the 90% confidence level. 
The corresponding upper limit product of branching fractions is B(B + —> 
© ++ p) x £>(© ++ —■> pK + ) < 9.1 x 10 -8 at the 90% confidence level. 

One theoretical conjecture [7] suggests that a possible glueball resonance Q 
(e.g., fj( 2220) [23]) with mass near 2.3 GeV/c 2 may contribute to the M pp 
threshold peaking behavior for the ppK + mode. Since the M pp mass resolution 
is about 10 MeV/c 2 , we scan through the 2.2 GeV/c 2 < M pp < 2.4 GeV/c 2 
mass region with a 20 MeV/c 2 wide window. The largest upper limit on the 
yield is found to be 18.9 at 2.21 GeV/c 2 . We use this data set to set an upper 
limit on the product of branching fractions of B(B + —» QK + ) x B(Q —> pp) < 
4.1 x ICG' at the 90% confidence level for a narrow glueball state with mass 
in the 2.2 - 2.4 GeV/c 2 range. The theoretical expectation is around 1 x 10 -6 . 

In summary, using 152 xl0 6 BB events, we measure the mass and the angu¬ 
lar distributions of the baryon-antibaryon pair system near threshold for the 
ppK + , ppKg and pKn~ baryonic B decay modes. A quark fragmentation in¬ 
terpretation is supported, while a gluonic resonant state picture is disfavored. 
Searches for a B meson decaying into the © + pentaquark or a glueball in the 
above related modes give null results. 
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Fig. 1. Differential branching fraction for (a) ppK + , (b) ppKg, and (c) pAn~ modes 
as a function of baryon-antibaryon pair mass. The shaded distribution shows the 
expectation from a phase-space MC simulation with area scaled to the signal yield. 
A charmonium veto has been applied in (a) and (b). 
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Fig. 2. Distributions of A E and M^ c , respectively, for (a) and (b) ppK + , (c) and 
(d) ppKg, and (e) and (f) pAir~ modes with baryon-antibaryon pair mass less than 
2.85 GeV/c 2 . The solid, dotted and dashed lines represent the combined fit result, 
fitted signal and fitted background, respectively. 
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Fig. 3. Branching fraction vs. cos 9 P in the baryon-antibaryon pair system for (a) 
ppK + , (b) ppKg, and (c) pAiT~ modes, (d) The proton angular distribution of the 
pir~ system against the A direction in the pAn~ mode, (e) ppK + background yield 
vs. cos 9 p . 



Fig. 4. Simple quark diagram for an illustration of the B + —>■ ppK + decay. 


13 











































Fig. 5. Distributions of M^ c and A E for the ppKg mode with 1.53 
GeV/c 2 < M pK o < 1.55 GeV/c 2 in the region where the 0 + pentaquark is ex¬ 
pected. The curves represent the fit projections. 
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